Metal-insulator transition in an aperiodic ladder network: an exact result 
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We show, in a completely analytical way, that a tight binding ladder network composed of atomic 
sites with on-site potentials distributed according to the quasiperiodic Aubry model can exhibit a 
metal-insulator transition at multiple values of the Fermi energy. For specific values of the first 
and second neighbor electron hopping, the result is obtained exactly. With a more general model, 
we calculate the two-terminal conductance numerically. The numerical results corroborate the 
analytical findings and yield a richer variety of spectrum showing multiple mobility edges. 

PACS numbers: 73.20.Jc, 72.15.Rn 



in 



CD 



O 

> 

o 

(N 

O 
00 
O 



X 



The classic problem of electron localization in low di- 
mensional quantum systems has remained alive over the 
last fifty years, since its proposition in 1958 by Ander- 
son [ll| . It is now well known that in one dimension even 
for arbitrarily weak disorder (almost )all the one-electron 
states are exponentially localized P, Q and one never en- 
counters mobility edges, that is, energy eigenvalues sep- 
arating localized (insulating) states from the extended 
(conducting) ones. Last couple of decades however, have 
witnessed an interesting twist in the canonical wisdom 
about localization problem through the advent of the so 
called correlated disordered models in one dimension. It 
is now established that even in a one dimensional chain 
of atomic sites, short range or long range [1, 0| po- 
sitional correlations between the constituents can induce 
unscattered (extended) wave functions at specific values 
of the electron-energy [1, |3l or even a crossover from lo- 
calized to extended states [sl, ^ . 

The search for a mobility edge in one dimensional lat- 
tices has always been an intriguing problem that started 
earlier with certain aperiodic lattices which bridge the 
gap between a perfectly periodic lattice and a completely 
random one. A famous example is the single band Aubry- 
Andre (AA) model [7| in which the on-site potential in a 
1-d chain of lattice constant a is described deterministi- 
cally by, e„ = A cos{Qna) with Q being an irrational mul- 
tiple of TT. It is to be appreciated that an AA-model repre- 
sents a certain class of almost periodic lattices which are 
much different from conventional randomly disordered 
systems, and displays a special kind of 'order', called the 
quasiperiodic order. The system however lacks transla- 
tional periodicity. In a one dimensional chain with near- 
est neighbor hopping integral t for example, and with 
such a potential (also called a Harper potential) the am- 
plitude /„ of the wave function at the nth site of the 
lattice can be obtained from the eigenvalue equation 



[E -\cos{Qna)] /„ = t(/«+i + /«-i) 



(1) 



In the AA-model, with X > 2t all the single particle 
states arc exponentially localized with Lyapunov expo- 
nent Zn(A/2t), while \ <2t makes the states (all of them) 
extended 0, Q , thus exhibiting a metal-insulator transi- 
tion in parameter space. No mobility edges exist for this 
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duality 
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Another premier feature of this model is its self 
which means that if one makes the transforma- 
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then the coefficients gm satisfy 



[E - 2tcos{Qma)] = -{gm+i + ffm-i) (3) 

The (7„j-equation is exactly of the same form as Eq. ([T| 
with the roles of t and A interchanged. It can be shown [8| 
that if the eigenstates given by gm are localized in recipro- 
cal space, then the eigenstates of Eq. ([1]) will be extended 
in real space, and vice versa. Incidentally, several vari- 
ants of the Aubry model have also been examined to de- 
tect the signature of mobility edges even in 1-d [i,[il[ii|. 

In this letter we investigate the electronic spectrum 
of an AA-ladder network built by fixing two identical 
AA-chains laterally (see Fig. [1]) . The motivation behind 
the present work is twofold. First, we wish to investi- 
gate if the interplay of the quasi-one dimensional struc- 
ture of the network and the AA duality, which is still 
preserved, leads to any possibility of a metal-insulator 
(MI) transition even within the standard form of the 
Harper potential. If it is true, then a ladder network 
such as this, could be used as a switching device, the 
design of which is of great concern in the current era 
of nanofabrication. Interestingly, research in AA-models 
in 1-d and its variations have been rekindled recently in 
the context of potential design of aperiodic optical lat- 
tices [H, [l^l • Therefore, the question of the existence of 
MI transition in a system with 'pure' AA (Harper) poten- 
tial can be addressed with a renewed interest. Secondly, 
the ladder networks have recently become extremely im- 
portant in the context of understanding the charge trans- 
port in double stranded DNA [3, ^M- Experimental 
results on DNA transport report wide range of behav- 
ior, from almost insulating 16], semiconducting (l7j to 
even metallic [isj . that can be attributed to many ex- 
perimental complications, such as the preparation of the 
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sample, sample-electrode contact etc. In addition to this, 
Mrevlishvili [l^ experimentally observed oscillations in 
the specific heat of DNA structures at low temperatures, 
results that have been subsequently explained by Mor- 
eira et al. [2^ considering a quasiperiodic sequence of 
the nucleotides. The results of reference [131, compare re- 
markably well with numerical results obtained for Ch22 
human chromosomes. It is to be appreciated that the 
helical structure of the double stranded DNA is expected 
to affect the periodicity of the effective site potentials on 
the ladder, and introduce incommensurate periods in the 
system. In view of this, the examination of the electronic 
spectrum of a ladder network comprising of aperiodically 
varying site potentials might throw new light into the be- 
havior of electrons, both in the context of basic physics 
and possible technological applications including DNA 
devices. 

We adopt a tight binding formalism, incorporate near- 
est and next nearest neighbor hopping inside a plaquette 
of the ladder and show that such a system exhibits a re- 



1 or 2. Using this basis, our task boils down to obtain 
solutions of the difference equation 

(£;i-e„)f„-t(f„+i+f„_i) (7) 

At first, we proceed to show the existence of multiple 
mobility edges in such a system in an analytically exact 
way. For this, we choose td — ti, and make the follow- 
ing transformation to the reciprocal space for each arm 
of the ladder and arrive at a difference equation in the 
reciprocal space, 

[{E ~ 2tiC0s{Qma)}I - {2t( cos((3ma) + 7} CTx] gmj 
= (A/2)I(g„+i.j+g™_ij) (8) 

where, CTx is the usual Pauli matrix. We now diagonalize 
the (Tx matrix by a similarity transformation using a ma- 
trix S, and define \(j)m) = S\gm) The difference equation 
([5]) now decouples into, 

(E + 7)0„,,2 = (A/2)((/)„,+ i.2 + (l)m-l,2) (9) 
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FIG. 1: Schematic view of a ladder attached to two electrodes. 

entrant MI transition. Most interestingly, for a selected 
set of the Hamiltonian parameters we are able to provide 
an exact analysis which shows the existence of mobility 
edges. We begin by referring to Fig. [TJ The Hamiltonian 
of the ladder network is given by. 



H = ^ CncJiCn + t ^ c]jC„+i + h.C. 



where 
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are the annihilation (creation) 
operator at the nth site of the jth ladder, e„_i = e„_2 = 
A cos{Qna) is the on-site potential at the nth site of the 
jth ladder, 7 is the vertical hopping between the nth sites 
of the two ladders, ti is the nearest neighbor hopping 
integral between the nth and the (n -t- l)th sites of every 
arm and td is the next nearest neighbor hopping within 
a plaquette of the ladder (see Fig. [1]). 

We describe the system in a basis defined by the vector 



fn,2 



(6) 



where, fnj is the amplitude of the wave function at the 
nth site of the jth arm of the ladder, j being equal to 



[E -J - Acos{Qma)] (j)m,i = {X/2){(j),n+i,i + 

(10) 

Here, 4>m,i and 4>m,2 are the elements of the column 
vector and A = 4i/. It is interesting to observe 
that, the Eq. ^ above corresponds to a perfectly or- 
dered chain with nearest neighbor hopping integral equal 
to A/2 in the reciprocal space. This implies that, for 
—A — 7 < E < A — 7 we have a gap less continuous 
spectrum in the reciprocal space. Eq. (jlOp on the other 
hand, represents the familiar single band Aubry model 
for which all states are localized or extended if A > A, 
or, A < A respectively. We can now extract information 
about the nature of eigenfunctions by considering the two 
Eqs. (O and pU]) simultaneously. 

Case I: \E + j\ < X and A > A 

We focus on the pair of Eqs. ([9]) and (fTO)) . When E 
lies within this range, we are within the 'continuous band' 
of extended states (in the reciprocal space). This means 
that the density of states corresponding to Eq. ^ is non- 
zero at all energies lying within this range and therefore 
4>m,2 7^ irrespective of the choice of A. Therefore, 

<S'2l5m,l + S'225m.2 7^ (11) 

for all m in dual space. This implies that, 

lim [^2l/„,l + S22fn,2] = (12) 
n — *oo 

in real space. Sij are the elements of the matrix S. If 
the average density of states corresponding to Eq. (fTO]) is 
non-zero, then Eq. (jlOp tells us that lim„i^oo (f'm,! — 0, 
as we have chosen A > A [8]. This, implies that, 

lim [Siigm,i + S'i25m,2] = (13) 

m — ^oo 

in dual space, and, 

Sllfn.l + Si2fn.2 + (14) 
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for all values of n in real space. Now, the density of 
states for an Aubry model is non-zero in the immedi- 
ate neighborhood of £' = 7 (the band-center) Q. We 
can therefore definitely say from Eq. (|12p and Eq. ([H]) 
that, both /„ 1 and /„_2 will be non-zero for any arbitrary 
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FIG. 2: g-E (red color) and p-i? (blue color) curves for a 
ladder of total number of rungs 60. (a) 7 = and (b) 7 = 3. 
Other parameters are, Q = {1 + \/E)/2, td = U = 3, eo = 0, 
to = 4 and A = 4. We have chosen c = e = h = 1. 

value of n. This ensures that all the states for the Aubry 
ladder will be extended at the center of the band. How- 
ever, interesting changes are observed as one looks away 
from the band-center. It is known [3, @| that the density 
of states of an infinite one band Aubry model is highly 
fragmented. This means, we shall have zero values of the 
density of states (i.e., no state at all) scattered through- 
out the spectrum. Whenever the density of states is zero, 
we shall have 



SllfnS + Si2fn,2 — 



(15) 

for all n. Considering Eq. I|12p and Eq. (|15p together, 
it becomes quite clear that, both and fn,2 have to 
be equal to zero as n oo in order that the Eqs. 
and are simultaneously satisfied. This means that 
the eigenstates are localized away from the center of the 
band. We thus have extended wave functions at the 
band-center flanked by the exponentially localized states 
on either side for —A — 7<i?<A — 7 and A > A. 

It is to be noted that in the true quasiperiodic limit 
the spectrum of an Aubry model exhibits more than one 
sub-bands separated by global gaps. Within each sub- 
band one has a highly fragmented band structure with 
infinitesimal energy gaps (for an infinite system). So, in 
principle following the argument given above, one should 
encounter an infinite number of mobility edges. How- 
ever, for realistic systems electron-electron interactions 



or the lead-sample connection will broaden the energy 
levels and the infinitesimal gaps will not persists. Only 
mobility edges which reside in the vicinity of the finite 
gaps separating the sub-bands will survive. 

Case II: \E + ^\> X and A > A 

In this energy regime E lies outside the band corre- 
sponding to the ordered system (Eq. (O), the correspond- 
ing density of states is zero (as there are no states at all). 
One then has. 



S2lfn,l + S22fn,2 — 



(16) 



for all n. If, on the other hand, the density of states 
corresponding to Eq. pop is non-zero, then 



S'll/n,! + Si2fn,2 7^ 



(17) 



for all n. Therefore, from Eq. (|16p and Eq. (|T7)) we ob- 
serve that, fn,i and fn,2 both remain non-zero for all 
values of n. This means, the eigenstates are extended for 
these values of energy. Thus, for A > A, the eigenvalue 
spectrum for the Aubry ladder, in real space, exhibits 
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FIG. 3: g-E (red color) and p-E (blue color) curves for a gen- 
eral AA-ladder of total number of rungs 60. The parameters 
are, 7 = 3, g = (1 + V5)/2, ta = 1, t; = 2, eo = 0, to = 4 and 
A = 4. We have chosen c = e = h = 1. 

the existence of localized and extended states separated 
by mobility edges and a re-entrant metal-insulator tran- 
sition is clearly visible. One can follow a similar chain of 
arguments to show that all states will be localized in real 
space for A < A. There are no mobility edges here. 

For a more general choice of the Hamiltonian param- 
eters (with ti ^ til E^nd e„.i 7^ £71,2)) analytical ap- 
proach becomes difficult. We have numerically calcu- 
lated the conductance of an Aubry ladder with various 
sets of parameters using a Green's function formalism. 
Though one has true quasiperiodicity only in an infinite 
system, for which we already have given an analytical 
proof for the existence of MI transition (for a special 
set of parameters), it is known that even finite systems 
grown following a quasiperiodic order are capable of ex- 
hibiting the localization effects [1, For example, 
in finite laboratory-grown Fibonacci multilayers experi- 
mental evidence of localization of light has already been 
reported [2l| . Therefore, though for an ideal infinite sys- 
tem several exotic spectral features may not be unlikely. 
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in our cases of interest, as we have worked in the param- 
eter regime where the single band AA-model has expo- 
nentially localized states only, expecting the localization 
fingerprints in a finite AA-ladder is quite legitimate. In 
the present calculation MI transition and mobility edges 
are found even when the simplification in the values of 
the hopping integrals are not made. 

To calculate the conductance, a finite Aubry ladder 
is attached to two semi-infinite one-dimensional metal- 
lic electrodes (Fig. [Ij, described by the standard tight- 
binding Haniiltonian and parametrized by constant on- 
site potential eo and nearest neighbor hopping integral to. 
For low bias voltage and temperature, the conductance g 
of the ladder is determined by the Landauer conductance 
formula [13] 9 = {2e^/h)T where the transmission prob- 
abihty T is given by [H T = Tr [TsGITdGI]. Ts and 
r D correspond to the imaginary parts of the self-energies 
due to coupling of the ladder with the two electrodes 
and Gl represents the Green's function of the ladder. 
All the informations about the ladder-to-electrodes cou- 
pling are included into these two self-energies as stated 
above and are described by the Newns- Anderson (NA) 
chemisorption theory [2^ . [24 | . For the sake of simplicity. 



here we have assumed that the entire voltage is dropped 
across the ladder-electrode interfaces and this assump- 
tion doesn't significantl y aff ect the qualitative aspects of 
the g-E characteristics [25j. In Fig. [5] we present the 
behavior of the conductance for the cases when 7 = 0, 
and ti = td, and the general case for a non-zero 7, and 
ti ^ td is shown in Fig. [3l In every case the pictures of 
the density of states are superposed to show clearly that 
we have eigenstates existing in energy regimes for which 
the conductance is zero. This spikes of localized eigen- 
states and the transition from the conducting (high g) to 
non-conducting phase is clearly visible in each case. 

Before we end, it should be pointed out that though 
the results presented in this communication are for zero 
temperature, they should be valid even for finite temper- 
atures 300 K) as the broadening of the energy levels 
of the ladder due to its coupling with the electrodes will 
be much larger than that of the thermal broadening [23] ■ 
The inter ladder hopping 7 will shift the spectra corre- 
sponding to Eqs. ^ and PU]) relative to each other, thus 
making it possible, in principle, to tune the positions of 
the mobility edges. This aspect may be utilized in de- 
signing a tailor made switching device. 
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